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Abstract. In this paper, we show that reflecting Brownian motion in any bounded domain 
D can be approximated, as fc — > oo, by simple random walks on "maximal connected" subsets 
of {2~^U^) n D whose filled-in interiors are inside D. 
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1. Introduction 

We proved in a recent article [1] that reflecting Brownian motion in a domain D can be 
approximated by a sequence of random walks on subsets of (2~'^Z'^) fl D. We chose A^s 
in a "natural" way, to be described in a moment. Our main theorem in [1] was limited to 
only some domains D ("extension domains"). We also provided a counterexample showing 
that random walks on A^'s do not converge to the reflecting Brownian motion in D for 
some domains D. In this paper, we will show in Theorems 13.61 and 14.21 that reflecting 
Brownian motion on any domain can be approximated by a sequence of discrete-time as well 
as continuous-time random walks if the state spaces Dk for the random walks are constructed 
in a different "natural" way. 

The sets were constructed in [1] as follows. First, we found the maximal connected set 
consisting of line segments contained in joining neighboring vertices in (2~^Z'^) flD. Then 
we let Ak to be the set of vertices in (2^'^Z'^) fl D at the ends of these line segments. It turns 
out that the "correct" way (employed in the present article) to construct the state space for 
the random walk is to start with the maximal connected set consisting of cubes contained 
in with edge length and vertices in {2~^'L'^) fl D. Then we let Dk to be the set of 
vertices in {2~^l/) fl D which belong to these cubes. Intuitively speaking, A^ may penetrate 
very thin crevices in D. The simple random walk on A^ may spend a non-negligible amount 
of time in such branches of A^ but reflecting Brownian motion spends very little time in sets 
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with very small volume. Replacing edges in the construction of A^s with cubes eliminates 
the mismatch between the shapes of D and the approximating discrete set. 

The technical essence of the paper is Theorem 12.11 which shows that, in a sense, the Dirich- 
let form for reflecting Brownian motion can be approximated from below by discrete Dirichlet 
forms. This theorem and the remaining part of the proof of the main result are challenging 
because "naive" discrete approximating schemes for the Dirichlet form of reflecting Brownian 
motion do not work — see Example 12.21 

In the rest of the Introduction, we will review some basic facts about reflecting Brownian 
motion in non-smooth domains and elaborate on some of the points mentioned above. 

Reflecting Brownian motion in a bounded domain D in Mf^ is a symmetric Markov process 
that behaves like Brownian motion inside D and is "pushed" back along the "inward normal" 
direction at the boundary dD of D. It is a prototype of diffusions with boundary condition 
and can be used to study heat equations with Neumann and Robin boundary condition. 
It is also widely used in modeling, for example, in physics, queuing theory and financial 
mathematics. Reflecting Brownian motion has been studied by various authors using various 
methods, see [H Ej and the references therein. When D is a bounded extension domain (see 
next paragraph for its definition), reflecting Brownian motion X can be constructed as a 
strong Markov process on D starting from every point in D except a polar set. Every bounded 
Lipschitz domain is an extension domain. When D is a general bounded domain, reflecting 
Brownian motion can still be constructed on D but typically it is no longer a strong Markov 
process. In a recent paper [1], we have developed three discrete approximation schemes for 
reflecting Brownian motion in bounded domains, providing effective ways to simulate the 
process in practice. The first two approximation schemes are discrete-time and continuous- 
time simple random walks on grids 2~^l/'PiD inside D. For these two approximation schemes, 
we need to assume that D is a bounded extension domain. A counter example is given in [1] 
showing that these approximation schemes do not work for some bounded domains. However, 
the third approximation scheme developed in [1], called myopic conditioning, works for any 
bounded domain D. Myopic conditioning generates a continuous-time and continuous-space 
process and, therefore, it is not suited for computer simulations. The purpose of this paper 
is to develop a discrete-time and continuous-time simple random walk approximations on 
grids inside D that work for every bounded domain D. 

We now give a precise description of reflecting Brownian motion on bounded domains. 
Let d > 1 and D he a, bounded domain in W^. The Sobolev space W^''^{D) of order 
(1,2) is the space of L^(D)-functions on D whose distributional derivative V/ is also 
L^(Z})-integrable. It is well known that W^'^{D) is a Hilbert space under norm ||/||i,2 := 

+ \\^f\\h{D^^ ^- We define on W^^^{D) a bilinear form 



£U.9) = \j^Vf{x)-Vg{x)dx for f,ge W^'\D). 
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It is known (see, e.g., [3]) that {£, W^'^{D)) is a Dirichlet form on L'^{D; dx). When C{D) fl 
is dense in both (C(D), || ■ \\^) and in {W^^'^{D), \\ ■ ||i,2), W^^'^{D)) is a regular 
Dirichlet form on L'^{D;m), where m is the Lebesgue measure on D extended to D by 
setting m{dD) = 0. In this case, there is a continuous conservative strong Markov process 
X on D associated with [S, W^''^{D)), starting from quasi every point from D. The process 
is called (normally) reflecting Brownian motion on D. It is known (see, e.g.. Theorems 1 
and 2 on pages 13-14 of [13]) that {S,W^''^{D)) is a regular Dirichlet form on L'^{D;m) 
if D is starshaped with respect to a point in D or if D has continuous boundary. Note 
that (£", W^^'^iW^) is a regular Dirichlet form on L'^iW^] dx). Hence (£", ly^'^(D)) is a regular 
Dirichlet form on L^{D;m) if D is an extension domain in the following sense: there is 
a linear continuous operator T : W^''^{D) — W^''^{M.'^) such that Tf = f a.e. on D for 
every / G W^''^{D). Recall that a domain D is called a locally uniform domain if there are 
6 G (0, oo] and C > such that for every x,y & D with |x — ?/| < 6, there is a rectifiable 
curve 7 in D connecting x and y with length(7) < C|x — y\ and moreover 

min{|a; — z\, \z — y\} < Cdist(2;, D^) for every 2 G 7. 

A domain is said to be a uniform domain if the above property holds with 6 = 00. The 
above definition is taken from Vaisala [I5], where various equivalent definitions are discussed. 
Uniform domain and locally uniform domain are also called (e, oo)-domain and {e, (5)-domain, 
respectively, in [12] . For example, the classical van Koch snowflake domain in the conformal 
mapping theory is a uniform domain in M^. Note that every bounded Lipschitz domain is 
uniform, and every non-tang entially accessible domain defined by Jerison and Kenig in [TT] 
is a uniform domain (see (3.4) of [H]), while every Lipschitz domain is an (e, 5)-domain. It 
is proved in [12] that every locally uniform domain is an extension domain. However for 
general domain (£^, W^^'^(Z)), £^) does not need to be regular on L^{D;dx). A unit disk 
in with a slit removed is such an example. See page 14 of [I3] for an example of D due 
to Kolsrud with dD = dD such that the Dirichlet form (£^, ly^'^(D), £^) is not regular on 
L'^{D; dx). Nevertheless, for any domain D C M.'^, one can always find a compact regularizing 
space D that contains D as a dense open subset such that {S, W^''^{D)) becomes a regular 
Dirichlet space on L'^{D] rh), where rh is the Lebesgue measure on D extended to D by setting 
m{D \ D) = 0; see [8] and [2]. Let X be the associated conservative strong Markov process 
on D, which can also be called reflecting Brownian motion on D. Let X be the projection 
of X onto D. Since for any given time t > 0, Fm{Xt G -D \ -D) = 0, under the normalized 
Lebesgue measure on D, X and X have the same finite dimensional distributions. 

A key technical element of this paper is to show that, for any bounded domain D in M'^, 
there exists a sequence {<^j,j > 1} of bounded smooth functions on D that is dense in 
the Sobolev space W^''^{D), separates points in D and satisfies the property (12. ip described 
below. We can deduce from its existence that there is a metric p on D ("refinement of 
the Euclidean metric") which induces the same Euclidean topology inside D and has the 



4 



KRZYSZTOF BURDZY AND ZHEN-QING CHEN 



property that reflecting Brownian motion on D can be lifted as a strong Markov process 
on the p-closure D oi D. This enables us to show that the random walk approximation on 
grids whose filled-in interiors are inside D works for reflecting Brownian motions on arbitrary 
bounded domains. In this paper, we also provide a proof that any weak limit of random 
walks on grids inside D is a stationary symmetric Markov process (see Theorem 13.31) . This is 
a key step in proving that the weak limit is indeed the stationary reflecting Brownian motion 
in D, using a Dirichlet form approach. This claim was made in [1] but regrettably no proof 
was given there. 

The rest of the paper is organized as follows. In Section [21 we establish a result (Theorem 
12. ip regarding the Sobolev space IV^'^(D) that will play an important role in this paper. 
Though the result is purely analytic, we employ some probabilistic techniques in its proof. 
The proof that reflecting Brownian motion in any bounded domain D can be approximated 
by discrete-time random walk on grids inside D is given in Section [3l The corresponding 
result for continuous-time random walk approximation is presented in Section HI 

2. Energy form estimates 

Let D C be a domain (connected open set) that has finite Lebesgue measure. Fix an 
arbitrarily small Ci G (0, 1) and a point Xq G D. For each integer /c, let Ak be the family of 
all closed d-dimensional cubes Q C D with edge length 2~'', such that 

(i) the vertices of Q belong to (2~^'Z)'^, 

(ii) the distance from Q to dD is greater than Ci2~^, 

(iii) there exists a sequence of cubes Qi,Q2, ■ ■ ■ ,Qm = Q, satisfying (i) and (ii), and such 
that xo G Qi, and Qj fl Qj+i is a (rf — l)-dimensional cube, for all j = 1, 2, . . . , m — 1. 

Since D has a finite volume, there is some ^Tj such that Ak = for every k < ko- Using 
scaling if necessary, we may and do assume that Ak = ^ for k < 0. Let Dk = Uqea. 

Theorem 2.1. Suppose that D G M."^ is domain with finite volume and C\ G (0, 1), Xq G 
D and Dk 's are defined as above. There exists a countable sequence of bounded functions 
{^j}j>i C l^i'2(D) n C^{D) such that 

(i) {^j}j>i is dense in W^''^{D), 

(ii) {'■Pj}j>i separates points in D, 

(iii) for each j > 1, 

(2.1) limsup 2^(2-'^) {^A^)-vMf<l \y^A^)?dx. 

Proof. Step 1. First note that the Sobolev space (ly^'^(il'), || • ||i,2) is separable. This can 
be seen as follows. Let Gi be the 1-resolvent for the Dirichlet form (£^, Vr^'^(D)); that is, Gi 
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is the linear operator from L'^{D;m) to W^'^{D) uniquely defined by 

Si{Gif,g)= [ f{x)g{x)dx for eyeij g e W''\D). 
Jd 

Here Si{u, v) := S{u, v) + u{x)v{x) dx. It follows that GiL'^{D] dx) is dense in the space 
II ■ II 1^2) and that 

£i(Gi/,Gi/)= / fGif{x)dx< [ fixfdx. 
Jd Jd 

Since L'^{D;dx) is separable, there is a sequence {fk,k > 1} of bounded functions that is 
dense in L'^{D;dx). Consequently, {77'^ := Gifk,k > 1} is a sequence of bounded functions 



1,2 



that is dense in {W^'^{D), 

Theorem 2 on page 251 of [7j implies that for every function rj'^, there exists a sequence of 
functions {r]j,j > 1} C W^''^{D) fl G°°{D) with the property that hmj_^oo ||?7^' - ^7^111,2 = 0. 
Moreover, the proof given in [7] shows that we can choose r]^ so that sup^jg^, |r7^(x)| < 
Ssup^g^ \r]^{x)\. 

Step 2. Constants ci, C2, . . . may change value from one "step" to another in this proof. 

We will use a regularized version of the distance function defined in [T^ Theorem 2, p. 
171]. That theorem implies that there exist < Ci, C2, C3, C4 < 00 such that for every integer 
j there is a C°° function dj : D ^ (0, 2~^] with the following properties. 

(2.2) ci(dist(x, dD) A 2"-') < dj{x) < C2(dist(x, dD) A 2^^), 



(2.3) sup |Vrfj(x)| < C3, 



(2.4) sup 

xeD 



dj{x) ^ ^ dj{x) 
OXi dxr 



< C4 for 1 < i,m < d. 



By dividing dj by an appropriate constant, we may and will assume from now on that 
(122])- (Q hold with C2 = C3 = 1. 

The existence of functions dj follows essentially from [HI Theorem 2, p. 171]. The only 
difference between our claim and that in [HI Theorem 2, p. 171] is that [H] is concerned 
with the condition 

cidist{x,dD) < dj{x) < C2dist{x,dD), 

in place of our condition (12. 2p . The method of proof given in [T3] applies to (12. 2p if we 
subdivide cubes constructed in [TH Sect. 1.2, p. 167] with edges longer than into cubes 
with edge length 2~^ . 

Step 3. Let : R'^ [0, 00) be a G°° "mollifier" with support in the ball 5(0, 1/2) such 
that Jg^Q -^^2) ip{x)dx = 1. For r > let ipri^) = r~'^ip{x/r) and note that sup^, ipr{x) = Cir~'^. 

The function ipl-di) := -^ipriv) is G°°. It is supported in _B(0,r/2) and satisfies the 
condition /^(or/2) '^riy)dy — 0- Let || ■ ||i denote the norm with respect to the Lebesgue 
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measure restricted to D. Note that W^j'^Wi = C2r~^ and \\ipr V 0||i = A 0||i = HV'rlli/S. 
Consider x E D and let 

'"IIV-i,.,vo|K ||>/.',,,.,A01|,- 

The functions a|^(-) and a~(-) are probabihty density functions. Let and A" be in- 
dependent M'^- valued random variables with densities a^(-) and a~{-), resp. Let dz 
denote the integral with respect to the length measure on the line segment joining and 
A~. Clearly, the measure E dz is supported on 5(0, dj{x)/2). We will now show that it 
has a density bounded above by Csdj{x)^~'^. In other words, for every set K C D, 

E / lK{z)dz < C3dj{xY~'^m{K n 5(0, d^{x)/2)). 
Ja+ 

Clearly the functions a^(-) and cl~{-) are bounded by := Cidj{x)~'^. Consider any 
z G B{x,dj{x)/2) and small 6 > 0. The probability that at least one of the random points 
or A~ belongs to B{z, 26) is less than c^S'^a^ < c^d'^dj^x)"''' . 

For /c > 2, the probability of the event Fk = {A+ E B{z,62'') \ 5(^,52'="^)} is bounded 
by C76'^2'^^ax. Let G be the intersection of B{x,dj{x)/2) and the smallest cone with vertex 
A^ containing B{z,6). The conditional probability, given F^, of the event {A~ G G} is 
bounded by times the volume of G, hence, it is bounded by Csaxdj{x){2~''dj{x)Y~^ = 
cgaxdj{xY2^^^~'^\ Multiplying the two estimates and summing over k > 2 such that 
B{z,62''~^) does not contain B{x,dj{x)/2) gives the bound 



C75'^2'^^axC^axdj{xY2^^^-'^^ < Cio6'^-^dj{ 

52^<2dj{x) 



xY-''. 



Adding to this quantity CQ5'^dj{xY'^ (the estimate representing the case when at least one 
of the random points ^4+ or A~ belongs to B{z,25)) gives a similar bound cii5'^~^dj{xY^'^ ■ 
The last quantity is an upper bound for the probability that the line segment joining A'^ 
and A~ intersects B{z,25). Since lB{z,2S){z)dz < 26 with probability 1, we obtain 

E / lBizM^)dz < cu5''-'dj{xY-''25 = cuS'^d^ixY"'. 

JA+ 

This estimate holds for all z G 5(x, dj{x)/2) and all small 5 > so the density of the measure 



E dz is bounded by cizdj{xY 

We will also need the following version of the above estimate. Let ipr{y) = -§^^r{y), 
bti-) = (V'i(.)(-)VO)/||^:;^(,)VO||i and6-(-) = -(<■(.)(• )AO)/||VJ(.)AO||i. The functions 
6^( ■ ) and b~{-) are probability density functions. Let 5+ and B~ be independent M'^- valued 
random variables with densities fe^ ( ■ ) and b~{-). The measure E dz has a density 
bounded above by Ci4^dj{xY~'^- In other words, for every set K C D,E j^^ lK{z)dz < 
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cudj{xY~'''m{K nB{0, dj{x)/2)). We omit the proof because it is analogous to the one given 
above. 

Step 4- Consider a function 77 G W^''^{D) fl C°°{D) and for integer j > 1 and x G -D, let 

(2.5) Vj{x)= ^pd,{x){y)vix -y)dy- 

J B(0,dj(x)/2) 

We will show that r]j G W^''^{D) fl C°°{D) and r]j — )■ 77 in ly^'^(D) as j — )■ cx). Since r] and 
are C°° functions, so is rjj. 
We have 

2 



(2.6) r/j(x) = ( / '^dj(x){y)v{x - y)dy 

'B{0,dj{x)/2) 



< / i^d,{x)iy) dy vi^-y) dy 

J B{0,dj{x)/2) J B(0,dj{x)/2) 

<c,{d,{x)-'fd,{xY [ Vix-yfdy 

J B{0,dj{x)/2) 



< C2dj{x) r]{x- yfdy. 

J B{0,dj{x)/2) 

Suppose that z G B{x,dj{x)/2). Then 

dist{z,dD) > dist{x,dD) — \x — z\ > dj{x) — \x — z\ > dj{x)/2. 

Hence, 

djiz) > C3(dist(2, dD) A 2~^) > C3{dj{x)/2 A 2'^) = c^dj{x)/2. 

Therefore, c^dj{z) > dj{x)/2 and lB{x,dj{x)/2){z) < lB{z,cid,iz)){x). Assuming again that 
z G B{x,dj{x)/2), we obtain 

dist(2, dD) < dist(x, dD) + \x - z\ < dist(x, dD) + dj{x)/2 
< dist(x, dD) + dist(x, dD)/2 < 2dist(x, dD). 

Hence, 

djix) > C5(dist(x, dD) A 2'^) > C5(dist(^, dD)/2 A 2'^) > c^{dj{z)/2 A 2"^) = C5dj{z)/2. 
This implies that 

(2.7) c?,(x)-^lB(.,d,(.)/2)(2) < C5-'^(rf,-(2)/2)-'^lB(.,c4d,(.))(a;). 

For later reference we derive an inequality slightly more general than what is needed in this 
step. For a set Q C il* let Q = Ux-gq B{x, dj{x)/2). We combine ( 12. 6p and ( 12. 7p to see that 

(2.8) / r]j{x)'^dx < C2 / dj{x)~'^ / r]{x - yfdydx 

Jq Jq J B{0,dj{x)/2) 
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C2 / dj{x) 



Q 



r]{zYdzdx 



B{x,dj{x)/2) 




C2 / Idjix) lB{x,dj{x)/2){z)ri{z) dzdx 



<C2 c^%d,{z)/2Y''lBi,,a,{z)){x)dxr^{zYdz 

JQ JQ 

= ce r]{zfdz. 



In particular, the inequality applies to Q = D = Q. Hence, 
(2.9) j r]j{xfdx <C6 [ 7]{zfdz. 

For any x & D, j and 1 < i < d, 



D 



(2.10) 



d 



dxi 



i J B{0,dj{x)/2) 



'ipd,ix)(,y)v{x - y)dy 



d_ 

dxi 



iB(o,djix)/2){y)^d,ix){y)v{x - y)dy 



d 



B{0,dj{x)/2) 



d 



d 



i^d,{x){y) ^^(^ -y)dy+ j — {'^B{o4,{x)/2){y)^di{x){y)) vi^ - y)dy 



dx 



i^d,{x){y) -^vix-y) dy + 

B{0,dj(x)/2) Cl^i JB{0,dj{x)/2) 



(^-^i^d,(x){y)^ Vix - y)dy 



<2l / 'ipd,{x){y) ^v{x - y) dy] +2 1 / ( w-'^d,{x){y)) v{x - y)dy 

' B{0,dj(x)/2) (y^i I \J B{0,dj{x)/2) \OXi J 



< 2 



B{0,dj{x)/2) 



'4^d,{x){yfdy 



B{0,dj{x)/2) 



-^Vix -y)] dy 



+ 2 



B{0,dj{x)/2) 



d 

g^'^d,{x){y) ] ri{x-y)dy 



< cjdjixy^ 

Jb 



B{0,dj{x)/2) 



_d_ 

dxi 



r]{x -y)] dy + 2 



B{0,dj{x)/2) 



d_ 

dxi 



iJd,[x){y) ]Ti{x-y)dy\ . 



Recall that il)[.{y) := -^ipriy) is a C°° function supported in B(0, 1/2) with Jg^^ ^^^-^ %jj'^(y)dy 
0. It follows from our construction of dj that Yli=i \ '£rdji^)\ — d- We have 



(2-11) E 



1=1 



d_ 

dx, 



■i^dj{x){y) 



E 



d_ 

dx, 



dj{x) 



\i^d,{x){y)\ < < Cgdjix) 



-d-1 
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Let ati-) = V0)/||^^^.(,) VOlli and «-(■) = " ) A 0)/||^^^,(,) A 0||i. The 

functions ■ ) and a~( ■ ) are probability density functions. Let and A~ be independent 
M''- valued random variables with densities a^{-) and a~{-). We have 



(2.12) 



B{0,dj{x)/2) 



d 



< ciodj{x) ^ \E{r]{x - At) - r]{x - A,))\ 



< Ciodj{xy^E \ri{x - At) - r]{x - A") I < cio(ij(x)~^E / |V?7(x - z)\dz, 



where the last integral is along a line segment from A^ to A^ . By Step 3, the measure 
E dz has a density that is bounded above by Ciidj{x)'^~'^ and vanishes outside of the ball 

5(0,rfj(x)/2). In other words, for every set C E Xf| lK{z)dz < Ciidj{x)^-^m{K n 
B{^,d.j{x)/2)). It follows that 

(2.13) 

d 



' B{0,dj{x)/2) 

< ciodj{x)^^ 
This implies that 



B{0,dj(x)/2) 



g^ i^d,{x)iy) ) ri{x-y)dy 



< ciodj{x)~^E / |V?7(x - z)\dz 



|V?7(x — 2;)|cii(ij(x)"'^ '^dz = Ci2d 



B{0,dj{x)/2) 



d 



|V?7(x — z)\dz. 



B(0,dj(x)/2) 



Q^ '^dj{x){y) ) v{x - y)dy | < | cudjix)' 



I V?7(x — z)\dz 



B{0,dj{x)/2) 



< c,,dj{x)-'%{x)'' 



B{0,dj{x)/2) 



|V?7(x - z)\'^dz = cisdjix)-'^ 



|V?7(x — z)\'^dz. 



B(0,dj(x)/2) 



We combine this estimate with fl2.10p to obtain 

2 



_d_ 

dx, 



■r]j{x^ 



<C7dj{x) [ ( ^ri{x - y)\ dy + cudj{x) f \Vri{x - z)\'^dz. 

JB{0,dJx)/2) \CXi J JB(0,di(x)/2) 



'B{0,dj{x)/2) 

Summing over i yields 



|V?7j(x)|^ < Ci5dj{x)~ 



|V?7(x — z)\'^dz. 



B{0,dj{x)/2) 



The same argument that leads from (12.61) to (12. 9p gives 



\Vr]j{x)f dx < / \'Vri{x)\'^dx. 



D 



D 



This formula and (12. 9 p show that r]j G W^''^{D). We have pointed out earlier in the proof 
that rij E C^{D). 
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Let Ks = {x e D : dist(x, D") > e}. We will show that r]j — )■ r/ in W^''^{D) as j oo. To 
see this, fix an arbitrarily small 6 > and find e > so small that 

{r]{xy + |V?7(x)p) dx < S. 

Note that the integral in the above formula is over the set A'g^, not K^. Since <Z D and 
r] is C°°, we have 



lim / ({r]j{x)-7]{x)y + {\Vr]Jx)\'^-Vr]{x)\'^))dx = 0, 

because the integrand converges to uniformly. It suffices to show that there exists a constant 
Ci7 < oo, not depending on 6 or e, such that for large j, we have 

(2.14) / {r]j{xy + \Vr]j{x)\^)dx <cn6. 

By (ESD applied to Q = K^, 

/ rij{x)'^dx < Cq rj{z)'^dz < cq6. 

A completely analogous calculation gives 

/ \\/rij{x)\'^dx < Cq6. 

The last two estimates yield f l2.14p and complete the proof of the claim that ry^ — )■ 77 in 
W^'\D) as J ^ 00. 

Step 5. Recall the constant ci G (0,1) and sets {Dk,k > 1} from the beginning of this 
section. For each integer k > 1, let Bk be the family of all closed d-dimensional cubes Q G D 
with edge length 2"^^, such that (i) the vertices of Q belong to (2~'^Z)'^, (ii) the distance from 
Q to dD is greater than Ci2~'^, and (iii) the interiors of Q and are disjoint. Let M.i = Bi 
and let JUk C Bk consist of those cubes in Bk that are not a subset of any cube in Bk-i for 
k>2. 

We recall that for a. set Q C D, Q = [J^^q B{x,dj{x)/2). We claim that there exists 
M < 00 such that every point x E D belongs to at most M distinct sets of the form Q where 
Q G IJ^ -Mk- This claim can be proved in a way that is totally analogous to the proof of [TU 
Prop. 3, p. 169] so we omit the proof. 

Step 6. We have shown in Step 1 that we can find a sequence of bounded functions {?7^}fc>i 
in iy^'^(D)nC°°(Z)) that is dense in W^''^{D). Let {f]^}j>i be a sequence constructed from t]^ 
as in (12.51) . The family {^7^ }fcj>i is dense in W^''^{D) and consists of bounded functions. 
Let us relabel the family {'U^}k,j>i as {'^j}j>i- We see that the family {'■Pj}j>i consists of 
bounded functions in W^''^{D) fl C°°{D) and satisfies part (i) of the theorem. 

By adding an appropriate sequence of functions in C^{D), if necessary, we can assume 
that condition (ii) is satisfied by {^Pj}j>i- 
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We will show that (12 .ip holds for (pj for each fixed j > 1. Some functions ^pj belong to 
C^{D). It is easy to see that (12. ip holds for such functions. Hence, we will assume that ipj 
belongs to the family {Vj}k,j>i- Then there exists a function cp in W^''^{D) fl C°°{D) such 
that ipj was constructed from ip as in (12. 5p . 

Fix an arbitrarily small e > and find an integer R so large that || V(/?1d^||l2(£)) < e. We 
increase i?, if necessary, so that dj{x) < for x G V^. Since (pj is C°°{D), we have by the 
mean value theorem, 

hm sup 2^(2-'^) ^ {ip,{x)-ip,{y)f< j |Vv.,f. 

It suffices to find Ci < oo independent of e and R and such that 
(2.15) limsup2'=(2-'i) ^ {^^{x)-ip^{y)f <cie. 



k—^oo 



Recall the notation from Step 5. Consider a large integer k, i < k and Q G Aie. Suppose 
that 

(2.16) Yl (ip.ix) - ipjiy))' = a. 

a;,j/G(2-'=Z)'*nQ,|x-j/|=2-'= 

Let = 2^^~^^'^ and let {Qi, Q2, ■ ■ ■ ,Qn} be the family of all cubes such that Q„ G Bk and 
Qn C Q. Let a.„ be the maximum of {<Pj{x) — (pj{y)Y over all pairs x,y E (2~'^Z)'^ fl Q.„ such 
that |a; — y| = 2"*^. By the mean value theorem, there is some z in the line segment joining 
X and y in Qn such that 

(2.17) \Vip,{z)\>al/'2^ 
It is easy to check that 

N 

(2.18) d2'^-^Yan>a. 

n=l 

Step 7. We have 

= TT" ( / ^'i.w(^) -y)dy+ f [-^i)d,(x){y) ) ¥^(3; - |/)rf?/ 

OXj \ Imd OXm ./TCPd \ OXm 



d 



-V'lz/) dy 
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B{x,dj{x)/2) 



d 



dXr 



+ 



B{0,dj{x)/2) 



d d 



i^d,{x){y) ](p{x-y)dy. 



We estimate the first of tlie last two integrals using (12. lip as follows 



Bix,dj(x)/2) 



d 



(2.20) < 



B{x,dj{x)/2) 



dx 

cidjix)-"^-^ 



iJdj{x){x-y) 
d 



d 



dXr. 



■^{y) dy 



dXr 



vi.y) 



dy < C2dj{x) 



-d-i 



\\/(p{z)\dz. 



B{x,dj{x)/2) 



To estimate the second integral, we apply the same method as in the derivation of (I2.13p . 
Recall that i/jKy) = ^tjjriy) is a C°° function supported in B{0, 1/2) with J^^^ ^j^^ ^lj'^.{y)dy = 



0. The function ip'^iy) = ■§piipr{y) is C°°. It is supported in 5(0,1/2) and satisfies the 
condition €{y)dy = 0. Note that = Cgr"!, ||^; V 0||i = \\^', A 0||i = Wr\\i/2, 



c^r-\ ||^;'VO||i = ||^:'AO||i = W;\\i/2. We have 



_d d_ 

j dx 77 



■^d,{x){y) 



d f f d 



dXi \\dxr, 
d d 



dj{x) ^rf^,(^)(?/) 



dXi dxr 



dj{x) 'ip'd^(x)iy) + 



d 



dXr 



■dj{x) 



^^.(x))<,.)(y). 



This implies that 
(2.21) 



B{0,dj{x)/2) \dXidXm 



d d 



'ipd,(x){y) ](p{x-y)dy 



< 



B{0,dj{x)/2) 



B{0,dj{x)/2) 



_d__d_ 

dx 2 dx 

d 



dj{x) 4^d^(x)iy) vix-y)dy 



dxr 



d 



dj{x) ] ( Q^dj{x) ) ipd^^^){y) ] ^{x - y)dy 



Recall that a+(-) = (^^,(.)(-) V Q)/\W,^^,) V 0||i and «"(■) = - (V^^^. (.)(■) A Q)/W,^^,^ A 
0||i. The functions a^(-) and a^(-) are probability density functions. Let and A~ 
be independent M'^- valued random variables with densities a^{-) and a~(-). Recall that 



I dxi 



dAx)\ < 1 and \-^-rr—dAx)\ < ddAx) ^ for all i.m and x. We obtain 



I dXi dXm ^ 

f 

B{0,dj{x)/2) 
< 



d d 



B(0,dj(x)/2) 
2 



dx, dxj'^'''^ J '^'^A-)^y^ ) '^(^ - ^)^^ 

c'rfj(x)-Vd^.(x.)(?/)) - y)dy 



< c,d,{xr' Eiifix - At) - ifix - A~)) 
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The same reasoning as in fl2.12l) -( l2T3l) yields 

d d 



(2.22) 



B(0,dj(x)/2) 



-d-1 



\Vip{x — z)\dz. 

B{0,dj{x)/2) 



We apply the same argument with ip" in place of ip'. Let 6+( ■ ) = ■ ^)/\\'^dj{x) ^^lli 

and 6~( ■ ) = — ■ ) A 0)/\\ip'^^f^^-^ A 0||i. The functions 6j( ■ ) and 6~( ■ ) are probability 
density functions that vanish outside the ball -8(0, dj{x)/2). Let and B~ be independent 
M*^- valued random variables with densities &J( ■ ) and b~{-). We have 



B{0,d,{x)/2) 



d 



d 



dxj'^"^^) [d^f'^"^^ ) <(x)(i/) ) - y)dy 



< cjd,{xr' \E{^{x - Bt) - ip{x - 5,7)) I 

< c^dj{x)-'^¥. \^{x - B+) - (f{x - B~)\ 



< C7dj{x)~^E / \V(p{x - z)\dz, 

where the last integral is along a line segment from B^ to B^ . By Step 3, the measure 
E dz has a density that is bounded above by csdj{xy~'^ and vanishes outside the ball 

B{0,dj{x)/2). In other words, for every set K C D, E lK{z)dz < Csdj{xY~'^m{K fl 
B{{),dj{x)/2)). It follows that 



(2.23) 



B(04j(x)/2) 



d 



dx„ 



Bx 



d 



djix) ) ( Q^dj{x) ) V^d^,(^)(y) ) <^{x - y)dy 



< Cjdjixy^E / |V^(x - z)\dz 



B2 



< Cjdj{x) 



Cgdj{x) 



B{0,dj{x)/2) 



|Vv?(x - z)\csdj{x) dz 



\Vip{x — z)\dz. 



B{0,dj{x)/2) 



We combine flCTj) . (KW^ . fl2:2T|) . (K22^ and fl2:23|) . and then we use H51der's inequality to 
see that, 



(2.24) 



d / d 



dxi V dx 



-VAX, 



< Ciodj{x) 



-d-l 



\V'^{z)\dz 



B(x4j{x)/2) 

< Cudjix)-'^-^dj{xy/^ 



1/2 



\V^{z)\^dz 



B(x,dj{x)/2) 
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= cnd.ix)-"/'-' ( [ \Vcp{z)\'dz 

\JB{x,dj{x)/2) 

Recall that Q = UxeQ -^(^' dj{x)/2). We will prove that 

(2.25) / |V(<^)P(x)dx>Ci22'=(2-d)^ 

Jq 

for some constant C12, where a is the constant defined in (I2.16p . If the inequality holds with 
C12 = 1 then we are done. So let us suppose that 

(2.26) / \V{ip)Wx)dx < 2'=(2-<i)a. 

Jq 

We combine this with 02.241) to see that for x G Q, 

d r d ^ ' " ^ 



dxi \ dx 



\UXm J \J B{x,dj(x)/2) 



'Q 

It follows from this and f l2.17p that the set oi x & Qn such that |Vv2j(a;)| > /2 contains 

a ball with radius greater than 

0130^^27 (rf,(x)"'^/^-i (2^(2-^) a)'/') = c,,al!\-"^2'''"^d,{xf^^\ 

and, therefore, it has a volume greater than 

[c,^a]l^a-^'^2^^/H,{xY'^+^Y = Ci,a^J^a-'"^2^^" 'H,{xY" . 

Hence 

|V<^,(x)|2dx> (0^2272) 'ci4a^/'a^'^/22^^'/2^,<x)'^'/2+d 
= Cr^al^'"^a-'"^2^^^+''''^^dj{xY"'^+\ 

and, therefore, 

(2.27) [ \V^,{x)\'dx = f2 [ \V^,{x)\'dx 

Jq n=i 

N 

> J2 c,,a\^'^'a-'/'2'^'+''/'^d,{xY'/'+'. 

n=l 

By the Holder inequality and (12.181) . 

n=l \n=l J 
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This and (12:271) give 

= cna2'^dj{xY{d,{x)2Y'^ > cisa2^''dj{xY > Ci9a2'^(2-d) 

It follows from this and (12. 8 P that, 
2'^'-"^ J2 (^.(^) - < C20 / \V^Ax))\'dx < C21 / \V^{x))\'dx. 

In view of f l2.25p - p.26p . we conclude that the last inequality is always valid. Recall the 
constant M from Step 5. Summing over all Q G IJ^<fe -^^5 Q obtain 

2'=(2-'^) J2 i^A^) - ^Mf < Mc,2 [ \Vip{x)\^dx < Mc^is. 

This shows that f l2.15p holds and completes the proof of the theorem. □ 



Example 2.2. Let Cl{D) be the family of bounded continuous functions on D with con- 
tinuous bounded first order derivatives. Using mean value theorem, it is easy to see that 
the inequality ( 12. ip holds for every ip e Cl{D) (in fact equality holds for such ip). However 
we will sketch an example, without proof, of a domain D such that Cl{D) is not dense in 
W^''^{D). The point of this example is to show that Theorem 12. II cannot be strengthened by 
adding an extra property that the functions {v^j}j>i belong to Cl{D). 
Let 





= {{Xi,X2) 


e ] 




-1 < Xi < 0, < X2 < 1}, 




= {iXl,X2) 


G ] 


: 


< a;i < 1, < X2 < 1}, 


Dn 


= {iXl,X2) 


G ] 


: 


— l/ra < xi < 1/n, 1/n < X2 < l/n + 5„}, 


dDt 


= {{Xl,X2) 


G ] 


: 


— l/n<xi < l/n,X2 = l/n + 6n}, 


dD- 


= {iXi,X2) 


G ] 


R^ : 


— 1/n < Xi < l/n,X2 = 1/n}, 


D 


= D_UD_ 


f u 


oo 

u 


oo 

D^\[j{dD:^udD;;)). 



n=2 n=2 



We choose 5„ > so small that -D„'s are disjoint. Consider a continuous function (p such 
that (p{x) = 1 for X G -D+ \ IJn>2 '^{x) = — 1 for x G -D+ \ IJn>2 and (p is linear in 
every The widths 6n of "channels" Dn can be chosen so small that ip G W^''^{D) and, 
moreover, iVy^P can be made arbitrarily small. 

We claim that the function p> described above cannot be approximated by functions t] G 
Cl{D) with arbitrary accuracy. The reason is that for any such 77, the oscillation of 77 in a 
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set Dn is arbitrarily small, for large n. Hence, in a neighborhood of (0,0), either \ip — iri\ is 
non-negligible on a non-negligible set, or \Wip — V77I is non- negligible on a non-negligible set. 
We leave the details to the reader because the claim made in this example is not needed for 
our main theorem. 

3. Invariance principle for reflecting random walk 

Let C be the algebra generated by functions {v^j}j>i from Theorem 12.11 over Q. By the 
same proof as that for Lemma 2.2 in [2], we have the following. 

Lemma 3.1. There exists a metric p on D which induces the same Euclidean topology 
inside D and such that the p-completion D of D is a regularizing space for Dirichlet form 
{S,W^-^{D)). Moreover, C is dense m Ch{D, \\ ■ ||oo)- 

Let m be the Lebesgue measure on D extended to D by setting m{D \ D) = 0. Then 
{S,W^''^{D)) is a strongly local regular Dirichlet form on L'^{D;m). Let X be the Hunt 
process on D associated with the regular Dirichlet form {S, W^''^{D)) on L^{D; m), which is 
continuous and has infinite lifetime. Denote by j the projection map from D to D. Then 
X := j{X) is a continuous Markov process taking values on D. In general X may not be a 
strong Markov process, as one can see from the example when D is the unit disk in with 
a slit (—1,0) X {0} removed. Both X and X can be called the reflecting Brownian motion 
on D. 

We will now discuss the relationship between reflecting Brownian motion X on D and a 
better known construction of reflecting Brownian on an arbitrary domain D. For an arbitrary 
bounded domain D in Mf^, Fukushima [S] used the Martin-Kuramochi compactification D* 
of D to construct a conservative continuous Hunt process taking values in D*. The process 
X* is associated with the regular Dirichlet form {S,W^''^{D)) on L'^{D*]m*), where m* is 
Lebesgue measure on D extended to D* by setting m*{D* \ D) = 0. Since each fi{x) = Xi is 
a function in W^''^{D), it admits a quasi-continuous extension to D*, which we still denote 
by fi. These functions induce a quasi-continuous projection map j* = (/i, ■ ■ ■ , fd) from D* 
to D. Then X' := j*{X*) is a continuous Markov process taking values on D, which is called 
reflecting Brownian motion on Z) in [2]. Both D and D* are regularizing spaces for the 
Dirichlet form {S, W^''^{D)) and so X and X' have the same finite-dimensional distributions 
under the initial distribution m. For x G -D, both X and X' starting from x behave like 
Brownian motion before they hit the boundary after a positive period of time. Consequently, 
X and X' have the same finite-dimensional distributions starting from any interior point in 
D. We can consider processes X and X' as maps from their underlying probability spaces 
into the space of continuous functions C([0, 00); M'^). Then the distributions of X and X' 
in C([0, 00); M*^) are identical either with initial distribution m or with initial starting point 
in D. In this sense, convergence of reflecting random walks to X or X' is an equally strong 
result. 
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Without loss of generahty, we assume that D contains the origin 0. Recall the definition of 
Dk from the previous section, connected. In this section, The graph with vertices {2~^'LY 
Dk and edges between nearest neighbors is connected. By abuse of notation, in this section 
we will use Dj, to denote {2~^ZY n D^. 

For X G -Dfc, we use Vk{x) to denote the degree of the vertex x in Dk- Let {X^^_2k,i = 
0, 1, ■ ■ ■ } be the simple random walk on Dk that jumps every 2"^^^ units of time. By definition, 
the random walk {X^^_2k,j = 0, 1,---} jumps to one of its nearest neighbors with equal 
probabilities. This discrete time Markov chain is symmetric with respect to measure rrik, 
where mk{x) = !^^il2~^'^ for x E Dk- Clearly rrik converge weakly to m on D. We now 
extend the time-parameter of |xj^2-2fc, J = 0, 1, ■ ■ ■ | to all non-negative reals using linear 
interpolation over the intervals ((j — 1)2~^^, j2~^'^') for j = 1, 2, ■ ■ ■ . We thus obtain a process 
= {Xf , t>0}. Its law with X^^ = x will be denoted by P^. 

For x,y E Dk, let x -H- y mean that x and y are at the distance 2~'^. Let Qk{x, dy) denote 
the one-step transition probability for the discrete time Markov chain |xj^2-2fe) J = 0, 1, ■ ■ ■ |; 
that is, for / > on D and x G Dk-, 



Qkf{x):= [ f{y)Qk{x,dy):=^ J2 fiv)- 

For / G C^(D), define 

Ckfix) := / ifiy) - f{x))Qk{x, dy) = ^ ^ {f{y) - /(x)), x G Dk- 

Then {X!:^_2^,) -^,^1^ Ckf{X^^_2^,), g^^_^,,j = 0,1,- ■ is a martingale for every / G 
C\D), where := a{X^, s<t). 

To study the weak limit of {X'',k > 1}, we introduce an auxiliary process defined 
by Y^^ := X^^2k^-^2-^'=' where [a] denotes the largest integer that is less than or equal to 
a. Note that Y^ is a time-inhomogeneous Markov process. For every fixed t > 0, its 
transition probability operator is symmetric with respect to the measure on Dk- Let 
J-'f := (t{Y^, s < t}. By abuse of notation, the law of Y'^ starting from x G Dk will also be 
denoted by Pj;. 

Note that Y^^ = X^ for every t of the form t = j2~^'^, where j is an integer. Moreover, 
supj>o \X^ — y^^l < 2~^. It follows that if the laws of one of the sequences {X^, > 1} or 
{y^, k > 1} converge to a limit on D([0, T],D) for some T, then the same holds for the other 
sequence. 

Theorem 3.2. Let D be a bounded domain in . Then the laws k > 1} of{Y^,k > 1} 
are tight in the space ©([0, T], D) for every T > 0. 
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Proof. Without loss of generality, we assume that T = 1. By P, Theorem 3.9.1] and 
Theorem 13.11 it suffices to show that for every g E C, {g{X^)^ is relatively compact in 
D([0, 1], M) with the initial distribution P^^. 

For each fixed k > 1, we may assume, without loss of generality, that Q is the canonical 
space ©([O, oo), D) and Y^^ is the coordinate map on Q. Given t > and a path u E Q, the 
time reversal operator is defined by 



(3.1) n{u){s) 



uj{{t-s)-), ifO<s<t, 
uj{0), iis>t. 



Here for r > 0, u{r—) := limg^r^is) is the left limit at r, and we use the convention that 
uj{0-) := uj{0). We note that 

(3.2) \imrtiuj){s) = u{t-) = n(w)(0) and limrt(a;)(s) = u{0) = n(a;)(t). 

Observe that for every integer T > 1, restricted to the time interval [0,T) is invariant 
under the time-reversal operator rj-. Note that 



[22'^t]-l 
i=0 



i=0 

is an {J-'t,t > 0}-martingale for every f EC. We have 

(3.3) f{Y,'^) - f{Y,') = ^Mf - i - Mf/„,^_) o for t E [0, T). 

For each M^'^ , there exists a continuous predictable quadratic variation process {M^'f)t. 
Note that (for example, see the page 214 of [9]) 

{M^^f), - {M^'f), = [Y. ^f^^u) - mrQkiY^', y)m,{y)du 

< 2d\\f\\Ut-s). 



Thus by Proposition VI.3.26 in [lO], {{M'''^)t}k>i is C-tight in D([0, 1],M). As ruk converges 
weakly to m on D, by [ini Theorem VI.4.13], the laws of {M'''f}k>i are tight in the sense 
of Skorokhod topology on ©([0,1], M) with the initial distribution P^fc- Since the laws of 
{M^-^, t E [0, 1], P^,Jfc>i are the same as the laws of {M^/_^^, t E [0, 1], P^Ja,.>i, it follows 
from (13.31) that {/(X^)}^^_^ and, consequently, {g{X'')^ is tight (and so relatively com- 
pact) in the sense of Skorokhod topology on ©([0, 1],M) with the initial distribution P^j.- ^ 



Let (X,P) be a subsequential limit of {(F^P^J; A; > 1} on B{[0,T],D). 
Theorem 3.3. (X,P) is a stationary symmetric Markov process. 
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Proof. Let (X,P) be a subsequential limit of F'^J; A; > 1} on D([0,T], 5), say, along 
a subsequence {nk,k > 1}. It suffices to show that the finite dimensional distributions of 
{X, P) are determined by a semigroup. Clearly, m is an invariant measure for X. For every 
t G [0, T], define a linear bounded operator on L'^{D; m) = L^{D; rh) by 

PJ = E[f{X,)\Xo], feL\D-m). 

We are going to show that {Ptit > 0} is a strongly continuous symmetric semigroup on 
L2(D;m). 

(i) We first show that each Pt is a bounded symmetric operator on L'^{D\m). For every 
f,g E Cb{D,p) and t > 0, it follows from the symmetry of {Y^,F^^) that 

fix)PMx)m{dx) = E \f{Xo)g{Xt)] = lim E^„^ [fiY,-')9iyn] 

£j L J fe— 5>oo '~ 

(3.4) = limE„„J^7(Fo"^)/(ll"^)]=EL(Xo)/TOl = / g{x)Ptf{x)m{dx). 
In particular, by taking = 1, we have 



(3.5) / Wtf\x)m{dx)= / f{x)m{dx) ioi f E CbiD, p). 

Jd Jd 

Note that Cc{D) C Cb{D, p) C L'^{D- m) and C^i^D) is dense in L'^{D] m). Hence ([33D holds 
for every / G L'^{D;m). Consequently, by the definition of Ft and Jensen's inequality for 
conditional expectation. 



(3.6) / {Ptfix)ym{dx) < / Pt(r)(x)m(dx) = / /(x)'m(dx). 

Jd 

Hence (13. 4p holds for every f,gE L'^{D, m); in other words, for each t > 0, is a symmetric 
contraction operator in L'^{D] m). 

(ii) Next we show that {Pt, t > 0} is a semigroup on L'^{D; m). For x = (xi, ■ ■ ■ , Xd) G D^, 
let Uk{x) := Y[f=i[xi — 2^'^~^ ,Xi + 2^^~^) be the half-closed, half-open cube centered at x. We 
define an extension operator : L'^{Dk,mk) — > L'^{D,m) as follows: for g G L'^{Dk,mk), 



(3.7) Ekgiz) :-- 



g{x) for z G Uk{x) with x G -Dfc, 

elsewhere. 



For f,g E Cc{D) and t of the form j2 by the uniform continuity, 

hm f f{x)E^^P:'g{x)m{dx) = \im E^^ [f{Y,')g {¥,')]= E^[fiXo)giXt)] 

fe— >00 j£) fe— >oo 

/(x)Pt5f(x)m((ix). 



D 



Note that 



(3.8) I {E^^Pt'"'gix)Ym{dx)= / (Pr^?(a;))^m„,(cia;) < / (?(x)^m„,(dx) 
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and that liiiifc^oo g{x)^mnf,{dx) = g{xYm{dx) for g e Cc{D). Since every / G 
L^{D; dx) can be approximated in L^-norm by a sequence {/fc, k > 1} C. Cc{D), we deduce 
from the last three displays and the Cauchy-Schwartz inequality that 



(3.9) 



lim / f{x)Enf^P^''g{x)m{dx) 
Jd 



f{x)Ptg{x)m{dx) 



D 



for every / G L'^{D; m) and g G Cc{D). 
We claim that for every t = j2~^', 



(3.10) 



lim 

k—^oa 



\En,P^'9{x) 



Ptg{x)\'^ m{dx) = for every g G Cc{D). 



D 



For any fixed x E D, there is r > so that B{x, 2r) C D. When k is large enough, there 
is a unique i/k G so that x G Ukiyk)- We denote this by 7rfc(x). Let q''(t,x,y) denote 
the transition density with respect to of simple random walk on killed upon leaving 
B{x,r). It follows from Donsker's invariance principle and the uniform Holder continuity 
[5], Proposition 4.1] for the parabolic functions of the simple random walk on 2'^%^ that 
g*^(t, 7rfc(x), nk{y)) converge locally uniformly in y G B{xo, r) to the transition density g(t, x, y) 
of Brownian motion X on R"^ with variance l/{2d) killed upon leaving B{x,r). For every 
e > 0, there is s > so that Jg^^ q{s, x, y)dy > 1 — (e/2). Hence for k sufficiently large, we 
have 

^^k{x)(Xs e dy) = q{s,x,y)mkidy) + Hk{dy), 

where Hk is a signed measure with |/ifc|(Dfc) < e. It follows from this and (13.91) that for 
9 e 



lim sup 
= lim sup 

fc— >oo 

< lim sup 

fc— >oo 

= lim sup 

fc— >oo 
= ^IIS'lloo- 



Pt'" 9{T^nk{x)) - / q{s,x,y)Pt-sg{y)dy 



D 



P;H^r.^7)K,(x))- / g(.,x,y)Pi_,5(?/)(i|/ 



D 



q{s,x,y)PiJ^,g{y)mk{dy) - j q{s,x,y)Pt^s9{y)dy 



q{s,x,y)En^P^J^^g{y)m{dy) - I q{s,x,y)Pt-sg{y)dy 
D Jd 



+ ^ 11^71 



Since e > is arbitrary, the above yields that {P^^'' 9{Tin^{x))\ > 1} is a Cauchy sequence 
and so it converges to some value u{x). This convergence holds for every x G -D, so we have 
by (13. 9p that u = Ptg a.e.; that is Enf.P^''g{x) = P^^ g{nnf.{x)) converges to Ptg{x) for a.e. 
X G -D. Hence by bounded convergence theorem fl3.10p holds for every g G Cc{D). 
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Abusing the notation a little bit, for / G L^(D;m), we define iikf as a function in 

L'^{Dk,mk) by 

(3.11) vr,/(a:) = ^^-(-);;7 . for x G D.. 

m(f/fc(a;)) 

Clearly vTfc o E"^ is an identity map on L'^^Di^^rrii^) and 

\'Kkf{x)\^mk{dx) < / |/(x)pm((ix). 
Since Cc{D) is dense in L'^{D]m), we have from (13.101) that 

(3.12) lim / |E.,P-.„,,(x) - f>.,(.)|^„(dx) = for every g e LHD;,n). 

It follows then for g G L^{D; m), 

Pt+s9 = lim E.^P^i^rvr^.^? = Jim (i?.,Pr^nJ(K,Prvr„J(7 = P^P.^?. 

fe— >oo fc— >oo 

This establishes the semigroup property of {Pt,t > 0}. 

We have now established that X is a stationary symmetric Markov process. □ 

The following result is needed in the proof of Theorem 13.61 

Lemma 3.4. In the above setting, for every f G C^{D), the process M/ := f{Xt) — 
/(Xq) — ^ ^f{Xs)ds is a ¥-square integrahle martingale. This in particular implies that 
{Xt^t < rD,P} is a Brownian motion killed upon leaving D, with initial distribution and 
infinitesimal generator ^A. 

Proof. Same as that for [H Lemma 2.2]. □ 

The following is Lemma 2.3 of pQ. 

Lemma 3.5. Let D be a bounded domain in and fix k > 1. Then for every j > 1 and 
feL\D,mk), 



if - Ql^f, f)L2{D,m^) < j{f - Qlf, f)L2{D,m^) < '^Hf " Qkf, f)L^D,, 



rrik)- 



We will say that "Zj is a Brownian motion running at speed 1/n" if Znt is the standard 
Brownian motion and we will apply the same phrase to other related process. 



By an argument similar to that in [H Section 2] but with Theorem l3.2l in place of [H Lemma 
2.2], and using Theorem 12.11 in the energy form argument in the proof of [H Theorem 2.4], 
we can establish the following theorem. 
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Theorem 3.6. Let D be a bounded domain in M'^ with m{dD) = 0. Then for every T > 0, 
the laws of {X'',F^^} converge weakly in C{[0,T], D) to a stationary reflecting Brownian 
motion on D running at speed \/d whose initial distribution is the Lebesgue measure in D. 
Consequently, for every T > 0, the laws of {X'',F^^} converge weakly in C{[0,T], D) to a 
stationary reflecting Brownian motion on D running at speed 1/d whose initial distribution 
is the Lebesgue measure in D. 



Proof. Fix T > 0. We know from Theorem 13.21 that the laws of (X , P^^) are tight in 
the space D([0,T],D). Let (X,P) be any of subsequential hmits, say, along (X'^J, Pm^,). By 
Theorem 13.31 and its proof, X is a time-homogeneous Markov process on D with transition 
semigroup {Ptit > 0} that is symmetric in L^{D,m). Let {Pt,t G 2~'^Z+} be defined by 
Pt^fix) := E^^[f{X^)]. For dyadic t > 0, say, t = and / G C, we have by Theorem O 

and Lemma [3.5[ 

t(/ - Ptf, /)L2m;m) = 7 (/ - Pt'f^ f)L^{D;nikA 

t ^ ' ' t J^oo J 



Jo ' ' 

< limsup ^jo2''^-''°(/ - Qk,f, f)L^D,m,.) 
j^oo Jo ^ 

= limsup2(2-'^)'=^i- 5^ {f{xf-f{x)f{y)) 

•' aiGDfc :y«x 



+ iimsup2(^-'^)^^i- ^ {fiyf-fi^)f(y)) 

y<^D^. x&Dkj-.y-^x 
= limsup2(2-'^)'^^i- Y U{x)-f{y)f 

x,yeDi,.:y^x 

Let {S,J^) be the Dirichlet form of X, or equivalently, of semigroup {Pt,t > 0}. That is, 

^UJ) = supi(/-PJ,/)^.(5^„) = limi(/-PJ,/)^,(5^^^ for/e^. 
Then for / G C, 

/) = sup \{f- PJ, <^l^ I V/(x)r dx. 
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This shows that / G J-". As C is dense in (W^''^{D), \\ ■ ||i_2) in view of Theorem I2.1[ we have 
W^'\D) C J" and 

SifJ) < for every / G W''\D). 

This, Lemma [3.41 and [H Theorem 1.1] (or [3l Theorem 6.6.9]) imply that J-" = W^'^{D) and 

S{fJ) = IW(a:)rrfx for / G ^-^'^(D). 

We deduce then that X is a stationary reflecting Brownian motion on D running at speed 1/d. 
This proves that converge weakly on C{[0,T], D) to the stationary reflecting Brownian 
motion on D running at speed 1/d. 

The last assertion comes from the fact that the projection map from {D, p) D is 
continuous. □ 

Remark 3.7. Note that for every x & D, X starting from x is a Brownian motion in D before 
hitting the boundary of D and the mass of X spreads immediately across the whole set D 
immediately after the clock starts, while X'^ starting from x runs like a simple random walk 
on 2~^Z before hitting the boundary of Dk- Using these properties, the weak convergence 
in Theorem 13.61 can be strengthened to show that (X'^,P^) converges weakly to (X,P^) for 
every interior starting point x E D. We leave the details to the reader. □ 

4. Continuous-time reflected random walk 

In this section, we show that reflected Brownian motion on D can be approximated by 
continuous-time random walks on grids. 

Let D he a. bounded domain in M"^ and be defined as in the beginning of Section [2l But 
in this section, X'' will be the continuous time simple random walk on D^, making jumps at 
the rate 2~^^. By definition, X'' jumps to one of its nearest neighbors with equal probabilities. 
This process is symmetric with respect to measure m^, where mk{x) = ^^^^^2~^'^ for x G D^. 
Note that converge weakly to the Lebesgue measure mon D, and recall that for x,y E D^, 
we write x ^ y ii x and y are at the distance 2~''. The Dirichlet form of X^ on L'^{Dk] rrik) 
is given by 

(4.1) s\f, /) = ^ 2-^^"'^'(/(^) - fiy))"- 

x,yeDk.x-ir^y 

Let Pj^^ denote the distribution of {X^,t > 0} with the initial distribution m^. 

Lemma 4.1. Assume that D is a bounded domain in M*^. For every T > 0, the laws of 
stationary random walks {X'^,P^^,/c > 1} are tight in the space 3{[0,T], D) equipped with 
the Skorokhod topology. 
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Proof. The proof is the same as that for [H Lemma 3.2] so we omit it. 



□ 



Theorem 4.2. Let D be a bounded domain in W^. Then for every T > 0, the stationary 
random walks on Dk converge weakly in the space D([0, T], D), as k ^ oo, to the sta- 
tionary reflected Brownian motion on D running at speed 1/d, whose initial distribution is 
the Lebesgue measure in D. 

Proof. Let P) be any of the subsequential hmits of {X^,F^J in the space B{[0,T],D), 
say, along X^^ . A similar argument as that for Theorem 13.31 shows that {Z,F) is a time- 
homogeneous Markov process and its transition semigroup {Pt,t > 0} is symmetric with 
respect to the measure m on D. Furthermore, by a similar argument as that in the proof of 
[H Lemma 2.2], the process Z killed upon leaving D is a killed Brownian motion in D with 
speed 1/d. Let {S,J^) be the Dirichlet form associated with Z, and let {Pt,t > 0} be the 
transition semigroup for X''. As X''^ converge weakly to Z in D([0, T], we have for every 
t>0 and every f,g e C(D), 

Let Ek : L'^{Dk;mk) — )■ L'^{D]m) and vr^ : L'^{D;m) — > L^{Dk,m) be the extension operator 
and restriction operator defined by (13. 7p and (13. lip , respectively. Then the last display can 
be restated as 

(4.2) hm (/, Ek^Pl''7rkjg)L^D;m) = (/, Ptg)L^D;m) for f,g e C{D). 

Note that \\Ekf\\L^D;m) = ||/||L2(Z)fe;mfe) for / G L'^{Dk\mk) and kg\\L^{Dk;mk) < hWh^iD^m) 

for g G L^{D] m). Since C{D) is dense in L^{D] m) and P^ and Pt are contraction operators 
on L'^{Dk;mk) and L'^{D;m), respectively, we deduce from (14. 2 p that 

(4.3) lim (/, Ek Pt'TTk g)L^D;m) = (/, Ptg)L^D;m) for every f,g E L'^{D; m). 
In particular, we have for every t > and f E C that 

(/, Ptf)L^D;m) = lim (/, Ek. P^' ITkJ) L^D;m) = lim (vTfc,/, Pt' TlkJ) L-^i^D^-mu)- 

Since C C Cb^D) and D is bounded, 

lim / \f{x) - 7rkf{x)\'^mk{dx) = for / G C. 
Hence we conclude that 

(4.4) (/, Ptf)mD;m) = lim (/, P''f)LHD,;m,) foT fEC. 

j-i>oo 

Thus we have for every t > and fEC, 

1 1 

t(/, / - Ptf)L^D;m) = ^m -(/, / - Pt'f)L^Dk;mk) 
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1 

-I 



< liminf sup -(/, / - Ps'f)LHD,:m,) 



j-s>oo Za ^ — ' 



< 

Thus 

£(fJ) = snplif-PtfJ)LHD;..)<^ [ \\/fix)\'m{dx) for every / G C. 

t>0 t ^CL Jd 

Since C is dense in the Sobolev space W^''^{D) with respect to norm || ■ ||i 2, it follows that 
J^dW^'\D) and 

£ifj) < ^J^\Vf{x)\Mdx) for every / G W'^'{D). 

Define 

S\f, 9) = j^j^ V/(x) ■ Vg{x)m{dx) for f,ge W''\D). 

Note that {S^, W^''^{D)) is the Dirichlet form for the reflected Brownian motion on D running 
at speed 1/d. On the other hand, as we have observed at the beginning of this proof, the 
process Z killed upon leaving D is a killed Brownian motion in D with speed 1/d. Therefore 
according to [U Theorem 1.1] (or [3i Theorem 6.6.9]), {S,J^) = {S^,W^''^{D)). In other 
words, we have shown that every subsequential limit of is reflected Brownian motion on 
D with initial distribution being the Lebesgue measure on D and with speed 1/d. This shows 
that X'' converges weakly on the space D([0, oo),D) to the stationary reflected Brownian 
motion X on D running at speed 1/d. □ 
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